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Dramatically  improving  the  performance  of  fuel  cell  systems  with  their  complex  heterogeneous  structures  involving 
electrocatalysts,  proton  conducting  membrane,  reactant,  and  interfaces  between  them  requires  understanding  the  fundamental 
chemical,  electrochemical,  and  physical  phenomena  at  the  heart  of  these  complex  materials  and  relating  these  fundamentals  to 
the  properties  and  performance  of  the  membrane -electrode  assembly.  Our  goal  is  to  develop  a  predictive  model  that  can  be  used 
to  estimate  the  changes  in  performance  upon  changes  in  the  design  and  which  can  be  used  to  monitor  performance  of  working 
fuel  cells.  Our  strategy  is  to  start  with  first  principles  quantum  mechanics  (QM)  and  to  develop  overlapping  simulation 
methodologies  in  which  QM  is  used  to  train  a  reactive  force  field  that  can  be  applied  for  large-scale  (millions  of  atom)  molecular 
dynamics  simulations  while  retaining  the  accuracy  of  QM.  The  results  of  molecular  dynamics  are  used  to  extract  a  coarse  grain 
or  mesoscale  description  useful  in  modeling  properties  at  much  larger  scales.  This  model  would  enable  the  conception, 
synthesis,  fabrication,  characterization,  and  development  of  advanced  materials  and  structures  for  fuel  cells  and  for  the 
associated  hydrocarbon  fuel  reformers  in  an  overall  fuel  cell  system.  We  illustrate  here  some  of  the  progress  toward  this  goal. 
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1.  Introduction 

Polymer  electrolyte  membrane  fuel  cell  (PEM-FC)  systems 
have  great  promise  for  providing  efficient  environmentally 
acceptable  portable  power  sources  for  transportation  and 
related  applications.  However,  despite  quite  significant 
investments  in  research  and  development,  current  systems 
lead  to  performance  inadequate  for  the  most  important 
applications  and  the  rate  of  progress  is  too  slow.  We  propose 
an  alternative  paradigm  for  attacking  this  problem  in  which 
theory  and  simulation  are  used  to  predict  from  first  principles 
the  structural,  transport,  and  electronic  properties  of  full 
PEM-FC  systems,  allowing  the  materials  and  configurations 
to  be  optimized  computationally  prior  to  experimental 
synthesis  and  characterization.  This  will  lead  to  current 
versus  applied  voltage  characteristics  for  various  materials 
choices,  the  proper  measure  of  overall  performance. 

Our  goal  is  to  determine  the  fundamental  processes  of 
PEM  fuel  cells  including: 

•  the  catalytic  reactions  at  the  anode, 

•  the  migration  of  the  protons  across  the  anode  -  electrolyte 
membrane  interface, 


•  the  transport  properties  of  these  protons  (and  fuel)  as  they 
traverse  the  electrolyte  membrane, 

•  the  migration  of  the  protons  onto  the  reactants  and 
intermediates  at  the  cathode,  and 

•  the  cathode  based  catalytic  reactions  to  produce  water. 

Such  studies  of  new  materials  and  configurations 
require  that  the  theory  and  simulation  be  based  on  first 
principles  (since  empirical  data  would  not  be  available), 
which  means  that  the  fundamental  description  of  the 
system  must  be  based  on  quantum  mechanics  (QM).  Until 
recently  it  would  have  been  completely  impractical  to 
consider  QM-based  studies  of  the  operation  of  fuel  cells. 
However,  a  recent  breakthrough  in  our  laboratory  provides 
the  promise  of  making  such  simulations  practical  and 
accurate. 

Recently,  we  have  developed  and  validated  the  ReaxFF 
reactive  force  field  to  describe  complex  reactions 
(including  catalysis)  [1,2]  nearly  as  accurately  as  QM 
but  at  costs  comparable  to  force  field  (FF)  based  molecular 
dynamics  (MD).  This  allows  us  to  consider  the 
fundamental  processes  for  realistic  descriptions  of  PEM 
fuel  cells. 
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These  ReaxFF  studies  can  be  carried  out  for  non¬ 
equilibrium  systems  in  which  current  and  protons  flow  so 
that  one  can  address  the  current -voltage  figure  of  merit 
relevant  to  FC  performance.  Each  of  these  processes  may 
require  simulation  models  on  the  scale  of  100  A,  requiring 
100,000-1,000,000  atoms,  but  this  is  now  feasible  with 
ReaxFF.  Such  ReaxFF  based  studies  of  the  dynamics 
would  enable  us  to  determine  the  fundamental  processes  at 
realistic  electrode -electrolyte  interfaces.  Currently  there 
is  essentially  no  information  about  these  interfacial 
processes,  which  we  believe  to  be  critical  in  determining 
fuel  cell  performance. 

We  have  already  validated  that  our  theory  and 
simulations  reproduce  some  of  the  known  properties  of 
PEM-FC  systems  by  simulating  well-studied  systems  (e.g. 
Nation  membrane  and  Pt  catalysts  for  anode  and  cathode) 
[3-5].  Now  we  are  ready  to  use  the  theory  and  simulation 
to  predict  the  properties  of  new  systems  in  which  catalysts 
are  designed  to  provide  lower  barriers  (overpotentials) 
with  reduced  effects  of  poisoning  and  in  which  the 
membranes  are  modified  to  optimize  proton  transport  at 
higher  temperatures  while  minimizing  transport  of  fuel, 
oxidant  and  products. 

This  approach  enables  a  new  strategy  for  fuel  cell 
development  in  which  we  consider  the  possibility  of 
combinatorial  computational  design  to  test  many  possible 
alloys  and  membrane  compositions  computationally  to 
find  the  best  predicted  performance  as  indicated  by  the 
overall  current -voltage  characteristics,  followed  by 
experimental  synthesis  and  characterization  on  these 
best  systems.  Elements  of  our  approach  include: 

•  Overlapping  simulation  paradigms  (figure  1),  in  which 
QM  data  are  used  to  train  first-principles  based  reactive 
force  fields  (ReaxFF)  that  are  employed  in  large  scale 
MD  simulations  to  follow  the  reactive  and  transport 
processes  at  the  heart  of  PEM-FC  systems. 

•  Developing  first-principles  (< ab  initio )  QM  for  proton 
transfer  studies  on  model  systems. 

•  Applying  newly  developed  density  functionals  QM 
(X3LYP)  for  estimating  the  thermodynamics  and 
kinetic  energy  barriers  of  proton  exchange. 

•  Determining  the  fundamental  chemical  reaction 
mechanisms  at  the  anode.  This  involves  understanding 
each  step  from  activation  of  the  fuel  to  forming  various 


Figure  1.  Overlapping  simulation  paradigms. 


intermediates  and  products.  We  will  determine  how 
this  depends  on  surface  defects  and  alloying.  In 
addition,  we  will  determine  how  each  step  is  affected 
by  solvent,  electric  field,  and  various  potential  poisons. 
Particularly  important  here  are  possible  poisoning  by 
heteroatom  components  in  the  fuel  (S  and  N)  and  by 
reaction  intermediates  (CO).  This  understanding 
should  allow  us  to  suggest  improvements  in  the 
anode  catalysts. 

•  Determining  the  fundamental  chemical  reaction 
mechanisms  at  the  cathode.  The  02  electroreduction 
at  the  cathode  is  kinetically  slow  (~  102  times  slower 
than  H2  oxidation  at  the  anode).  Typically,  the  platinum 
catalysts  consist  of  nanoparticles  supported  on  a  carbon 
electrode.  We  have  already  made  some  progress  in 
determining  the  detailed  mechanism  of  the  cathode 
reactions  on  the  Pt(lll)  surface  and  on  multi- metallic 
(Pt/Ni,  Pt/Co,  Pt/Mo,  etc.)  alloys.  We  plan  to  extend 
these  studies  towards  understanding  the  mechanism  of 
the  cathode  reactions  (including  barriers)  while 
including  the  effects  of  the  carbon  support  and  the 
presence  of  the  electrolyte.  We  expect  that  such  studies 
will  lead  to  suggested  improvements  in  the  structure 
and  composition  of  the  cathode  catalyst  materials. 

•  Developing  ReaxFF  and  hybrid  classical/ReaxFF 
molecular  dynamics  (hybrid  ReaxFF  MD)  to  model 
the  fundamental  fuel  cell  chemistry  while  keeping  the 
expense  of  the  calculations  relatively  low. 

•  Developing  MD  and  mesoscale  methods  for  determin¬ 
ing  solubility,  diffusion,  and  permeability  of  protons 
and  water. 

•  Determining  the  kinetics  for  transport  of  protons,  H20, 
fuel,  fuel  products,  and  02  through  the  PEM  and  how 
this  depends  on  the  nanostructure,  concentration  of 
water  and  salt,  pH,  temperature,  and  nature  of  the  PEM. 

•  Characterization  of  how  the  nanostructure,  heterogen¬ 
eity,  and  properties  of  the  PEM  depend  on  polymer 
architecture. 

•  Detailed  study  of  the  interfaces  between  the  electrode 
catalysts  and  the  PEM  with  an  emphasis  on  how  this 
affects  the  transfer  of  the  protons  between  the  PEM  and 
the  appropriate  reaction  intermediates  on  the  catalysts. 
This  is  likely  the  source  of  both  overpotential  and 
retarded  voltage  at  high  currents.  We  will  then  modify 
this  interface  by  changing  the  catalyst  and  membrane 
to  determine  new  materials  that  improve  catalyst 
performance. 

•  Developing  fuel  cell  models  that  allow  us  to  suggest 
new  leads  for  new  electrolytic  membranes  to  improve 
performance  or  to  change  the  operating  conditions  (e.g. 
temperature).  We  consider  new  alternative  membranes 
(such  as  the  imidazole  and  dendrimer  based  modified 
electrolytes  already  suggested).  Our  goals  are  high- 
proton  conductivity  at  appropriate  temperature,  low 
fuel  and  oxidant  crossover,  and  excellent  chemical  and 
mechanical  stability. 

•  Developing  a  full  device  fuel  cell  model  based  on  the 
understanding  of  the  fundamental  processes  (electrode 
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reaction  kinetic,  internal  transport  coefficients  for 
protons,  water,  and  fuel)  to  be  obtained  in  the  above 
studies.  We  focus  on  simulating  the  phase  separation, 
morphology,  interface  structures  of  composites  affect¬ 
ing  the  overall  fuel  cell  performances. 

For  the  most  rapid  progress,  it  is  essential  to  couple 
these  computational  efforts  to  experimental  tests  of  the 
predictions  with  the  results  of  subsequent  experiments  fed 
back  into  the  theory,  perhaps  leading  to  refined 
simulations  aimed  at  resolving  discrepancies  or  consider¬ 
ing  new  phenomena.  Ultimately  this  process  of  developing 
a  complete  computational  model  whose  predictions  have 
been  validated  by  experiment  will  allow  an  accurate 
engineering  model  to  be  developed  based  on  the  theory 
and  simulations  but  incorporating  all  the  results  obtained 
from  experiments.  Only  with  such  wide-ranging  collab¬ 
orations  between  experiment  and  theory  it  will  be  possible 
to  develop  fully  validated  theoretical  models  for 
developing  new  materials  and  fully  validated  engineering 
models  for  optimizing  performance. 

2.  Fundamental  theory 

In  order  to  predict  from  theory,  the  properties  of  new 
materials  (prior  to  synthesis)  first-principles-based 
methods  must  be  used.  However,  these  models  must  be 
adequate  to  predict  the  large  spatial  and  time  scales 
relevant  to  devices.  To  accomplish  this  we  use  an 
overlapping  hierarchy  of  computational  methods  involv¬ 
ing  different  scales  but  each  building  upon  the  elements  of 
a  more  fundamental  description.  Thus,  the  QM  is  used  to 
determine  the  mechanisms  of  catalytic  reactions  and  to 
predict  the  characteristics  observable  experimentally  (e.g. 
XANES,  EXAFS,  NMR).  In  addition,  the  QM  results  are 
used  to  develop  FFs  for  large-scale  MD  calculations. 
Similarly,  MD  results  are  used  to  extract  properties  for 
larger  length  and  time  scales  (e.g.  diffusion  of  H20  and 
protons). 

Additional  developments  are  required  in  the  QM  study 
of  the  role  of  interfacial  structures  and  applied  voltages  in 
electrocatalysis  and  the  use  of  ReaxFF  to  describe  the  flow 
of  protons  between  catalyst  electrodes  and  electrolyte  and 
within  the  electrolyte. 

2.1  Density  functional  theory  methods 

We  have  used  many  QM  methods  to  probe  chemical 
reaction  mechanisms  and  find  that  the  B3LYP  and  X3LYP 
[6]  flavors  of  DFT  generally  provide  adequate  accuracy  at 
modest  cost.  We  use  this  approach  to  characterize  the 
catalytic  processes  at  the  anode  and  cathode.  Our  strategy 
here  is  to  carry  out  QM  calculations  with  the  catalyst 
represented  by  clusters  of  —  35  atoms  [7]  to  determine  the 
mechanism  of  the  dispersed  catalyst  and  then  to  do  QM 
calculations  on  periodic  slabs  to  describe  (ordered) 
overlayers  on  infinite  surfaces.  For  the  finite  clusters,  we 
use  Poisson- Boltzmann  continuum  solvent  methods  to 


provide  an  accurate  description  of  solvent  polarization. 
We  have  applied  this  successfully  to  reactions  in  water  and 
sulfuric  acid,  where  it  leads  to  solvation  energies  (for 
neutral  species)  good  to  —0.5  kcal/mol  and  to  p Ka  values 
good  to  —0.5.  We  also  use  these  methods  to  describe  the 
elementary  steps  for  proton  transfer  through  the 
electrolyte. 

For  periodic  calculations,  we  prefer  the  SeqQuest  code 
being  jointly  developed  at  Sandia  National  Labs  (Peter 
Schultz)  [8]  and  at  Caltech,  and  the  Crystal2003  software 
[9].  Both  allow  the  use  of  Gaussian  basis  functions  and 
hard  or  soft  pseudopotentials. 

2.2  First  principles-based  reactive  force  field  ( ReaxFF ) 

The  QM  calculations  described  above  are  useful  for 
characterizing  idealized  electrochemical  processes  (invol¬ 
ving  simple,  ordered  structures,  zero  temperature,  etc.). 
However,  these  methods  are  too  slow  for  routine  studies  of 
the  fundamental  unit  processes  that  characterize  FC 
performance  under  actual  operating  environments  (i.e. 
realistic  electrode -electrolyte  interfaces,  finite  tempera¬ 
ture,  disordered  structures,  longer  time-scales).  MD 
calculations  using  FFs  provide  a  computationally  efficient 
way  of  describing  the  dynamics  with  large-scale 
simulations.  However,  FFs  have  not  generally  been 
capable  of  accurately  describing  chemical  reactions. 
Recently,  we  have  made  a  breakthrough  in  this  area  with 
the  development  of  the  ReaxFF  capable  of  reproducing 
high-level  QM  data  for  reactive  and  non-reactive  systems. 
ReaxFF  descriptions  have  already  been  reported  for  a  wide 
range  of  systems,  including  hydrocarbons  [10],  nitramines 
[11],  peroxides  [12],  aluminum  metal  and  aluminum 
oxides  [13],  silicon  metal/silicon  oxides  [14]  and  silicon 
carbides  [1],  copper/nickel/cobalt  interactions  with  carbon 
[2],  magnesium  and  magnesium  hydrides  [15],  lithium/- 
lithium  carbides  [16]  and  all-carbon  materials  [17]  while 
ongoing  developments  are  extending  the  methodology 
further  across  the  periodic  table  (figure  2),  including 
reactions  on  Pt,  Pt3Co,  Pt3Ni,  PtRu,  and  BiTeOv  surfaces. 
In  addition  a  similar  methodology  provides  an  excellent 
description  of  a  BaTi03  ferroelectric  materials  [18]. 

This  progress  demonstrates  that  a  wide  range  of 
reactions  and  reactive  systems  can  be  described  using 
essentially  the  same  FF.  ReaxFF  has  also  been  re-coded 
for  a  parallel  environment,  and  fully  reactive  simulations 
of  around  half  a  billion  atoms  have  been  reported  [19], 
suggesting  that  realistic  simulations  covering  the  full 
complexity  of  the  fuel  cell  environment  are  now  feasible. 

The  ReaxFF-description  involves  the  following  three 
types  of  interactions: 

a.  Environment  dependent  charge  distributions  on 
atoms.  The  charge  on  each  atom  is  described  with  a 
Gaussian  shape  to  provide  for  proper  shielding  as  the 
atoms  come  close  together.  This  allows  the  electro¬ 
static  interactions  of  bonded  atoms  to  be  included, 
eschewing  the  common  approach  of  “excluding” 
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Figure  2.  Atoms  currently  described  by  the  ReaxFF  method. 


bonded  atoms  from  the  electrostatics.  The  charge  of 
the  atom  is  allowed  to  transfer  to  other  atoms  in 
response  to  the  local  environment  of  each  atom  [as  in 
the  charge  equilibration  (QEq)  model].  Such  proper 
description  of  electrostatics  is  essential  in  describing 
systems  undergoing  oxidation  and  reduction,  as  in  fuel 
cells. 

b.  Bond-order  dependent  valence  terms.  All  valence 
terms  (bonds,  angles,  torsions)  go  to  zero,  as  the  bonds 
are  broken,  allowing  a  smooth  description  of  chemical 
reactions.  The  bond  orders  are  determined  only  by  the 
interatomic  distances  and  are  updated  every  iteration, 
allowing  ReaxFF  to  automatically  recognize  and 
handle  changes  in  connectivity  as  reactions  proceed. 

c.  General  non-bond  van  der  Waals  interactions.  To 
account  for  short  range  Pauli  repulsion  and  longer 
range  dispersion  interactions  we  use  a  universal 
vdW  term  that  applies  to  every  atom  pair, 
independent  of  bonding  (again  no  exclusions). 
These  terms  are  also  shielded  to  properly  describe 
the  changes  at  short  distances. 

The  following  guiding  principles  were  used  in  the 

ReaxFF  development: 

d.  Every  parameter  is  derived  solely  from  QM.  A 
critical  aspect  of  the  ReaxFF  is  that  it  is  based 
solely  on  ab  initio  QM  calculations  and  parameter¬ 
ized  to  reproduce  a  large  number  of  reactions 
designed  to  characterize  the  atomic  interactions 
under  all  possible  environments  an  atom  might 
encounter. 

e.  Transferable  potential.  Each  element  is  described 
by  just  one  atom  type.  That  is  there  is  one  set  of 
carbon  FF  parameters,  which  automatically  describe 
correctly  C  in  diamond  and  graphite,  organics  like 
ketones  and  allyl  radicals,  the  carbide  on  top  of  a 
Ni  surface.  This  is  necessary  since  the  reactions 
may  interconvert  between  all  of  these  forms. 

f.  No  predefined  reactive  sites.  Since  the  concept  of 
bonds  and  reactions  arises  naturally  in  terms  of  the 
interactions  (just  as  in  QM)  the  user  does  not 


predefine  where  and  when  reactions  are  expected  to 
occur.  Thus  ReaxFF  allows  unbiased  simulations  of 
reactive  systems. 

g.  No  discontinuities  in  energy  or  forces.  Even  during 
reactions,  ReaxFF  provides  a  continuous  energy  and 
force  description,  thus  allowing  proper  reactive  MD 
simulations. 

Recently,  we  have  developed  ReaxFF  potentials  cover¬ 
ing  the  following  aspects  of  PEM-FC  chemistry: 

(1)  Reactions  at  metal  surfaces.  To  enable  a  large-scale 
dynamical  description  of  the  chemical  events  at  fuel 
cell  metal  anode  and  cathode  we  have  developed 
ReaxFF  for  Pt/C/H/O,  Ru/H/O  and  Ni/C/H/O 
interactions. 

(2)  Oxidation/reduction  reactions.  To  ensure  that 
ReaxFF  describes  redox  reactions  accurately,  we 
have  fitted  ReaxFF  against  a  wide  range  of  metal 
oxidation  states. 

(3)  Hydrogen  transport.  The  approach  that  we  have  used 
to  predict  proton  transport  through  complex  mem¬ 
branes  such  as  Nafion  is  to  use  QM  to  determine  the 
barriers  for  H  migration  as  a  function  of  the  donor - 
acceptor  separation  and  then  to  combine  probabilities 
based  on  these  results  with  MD  studies,  allowing 
proton  transport  to  occur  when  appropriate  configur¬ 
ations  are  encountered. 

With  ReaxFF  we  do  not  need  to  separate  these 
processes,  since  the  reactions  can  occur  during  the 
dynamics.  To  validate  the  application  of  ReaxFF  to 
hydrogen  transport  through  electrolytes  figure  3  compares 
QM  and  ReaxFF  results  for  H-migration  between  an 
imidazole  and  an  imidazolium  cation  as  a  function  of 
nitrogen -nitrogen  distance  [20].  This  shows  that  ReaxFF 
properly  describes  the  distance  dependent  hydrogen 
migration  barriers,  indicating  that  ReaxFF  can  be  used 
to  describe  reactive  transport  through  fuel  cell 
electrolytes. 
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QM  °  ReaxFF 


Figure  3.  QM  and  ReaxFF  energies  for  hydrogen  migration  barriers  between  an  imidazolium  cation/imidazole  pair  as  a  function  of  nitrogen -nitrogen 
distance. 


Since  the  ReaxFF  applies  to  all  combinations  of  the 
elements,  the  ReaxFF  descriptions  of  separate  components 
of  PEM-FCs  can  be  used  on  the  composite  systems, 
allowing  the  ReaxFF  based  simulations  for  studying  the 
critical  chemical  events  at  PEM-FC  interfaces. 

2.3  Large-scale  molecular  dynamics  simulations 

MD  simulations  using  the  accurate  ReaxFF  will  provide  a 
first-principles  characterization  of  the  fundamental 
processes  (structures,  mechanisms,  energetics,  and  rates) 
and  mechanisms  in  a  FC  including: 

1.  Poly  electrolyte  membranes 

a.  Microstructure  of  hydrated  polymers  (Nafion®) 
and  new  alternative  membranes; 

b.  Proton  transport  in  hydrated  polymer  (Nafion®) 
and  new  alternative  membranes: 

(1)  Mechanisms; 

(2)  Role  of  water/water  transport. 

2.  Anode/electrolyte  interfaces 

a.  Obtain  atomistic  structures  of  the  interface 
between  the  anode,  electrolyte,  gas  (fuel); 

b.  Characterize  the  chemistry  at  the  three-phase 
interface  under  realistic  operating  conditions  (this 
includes  the  migration  of  the  proton  to  the 
electrolyte  and  the  detailed  reaction  steps  on  the 
catalyst). 

3.  Cathode/electrolyte  interfaces 

a.  Obtain  atomistic  structures  of  the  interface 
between  the  cathode,  electrolyte,  gas  (02); 

b.  Characterize  the  chemistry  at  the  three-phase 
interface  under  realistic  operating  conditions  (this 
includes  the  migration  of  the  proton  from  the 
electrolyte  and  the  detailed  reaction  steps  on  the 
catalyst). 

Many  of  these  applications  require  large-scale  MD 
simulations  on  10,000-1,000,000  atoms.  These  will  be 
based  on  our  recently  developed  computational  materials 
design  facility  (CMDF)  that  allows  a  mixture  of  paradigms 


and  scales  to  be  utilized  simultaneously  to  describe  the 
time  evolution  of  large  scale  complex  reacting  systems. 
This  framework  allows  the  facile  mixture  of  complex 
independent  software  codes  development  for  a  various 
paradigms  (QM,  ReaxFF,  classical  FF,  mesoscale  FF,  finite 
element  analysis)  to  be  used  independently  but  simul¬ 
taneously  on  the  same  problem. 

3.  Computational  modeling  of  electrode/electrolyte 
interfaces 

To  develop  the  capability  for  modeling  the  electrode/elec¬ 
trolyte  interface  in  a  PEM-FC  we  examine  the  properties 
of  Nafion  hydrated  polymer  as  membrane  and  Pt  catalysts 
as  the  electrode  material.  There  is  a  great  deal  of  data  on 
Nafion  and  Pt  for  both  cathode  and  anode  in  a  PEM-FC, 
but  little  microscopic  understanding  of  the  properties.  This 
availability  of  experimental  data  allows  us  to  check 
reliability  of  our  first-principles  based  model  by  compar¬ 
ing  to  experiment. 

3.1  Nanostructure  and  transport  properties  of  Nafion 

Although  great  progress  has  been  made  since  the 
pioneering  work  by  Gierke  on  the  Nafion  structure  in 
1980s  [21],  the  details  of  the  microstructure  of  hydrated 
Nafion  remain  highly  disputed.  The  main  function  of  the 
electrolyte  membrane  in  a  fuel  cell  is  to  provide  high 
proton  conductivity,  electrochemical  and  mechanical 
stability,  and  a  barrier  for  transport  of  electrons,  fuel, 
and  oxidant.  Despite  its  high  cost  and  restriction  to  work  at 
temperature  not  higher  than  80°C,  Nafion  remains  the 
material  of  choice  for  fuel  cell  applications. 

Nafion  is  a  polyelectrolyte  consisting  of  non-polar 
tetrafluoroethylene  (PTFE)  segments,  N  =  (CF2-CF2) 
and  polar  perfluoro sulfonic  vinyl  ether  (PSVE)  segments, 
P  =  (CF2  -  CF(0  -  CF2  -  CF(CF3))  -  CF2  -  CF2  -  S03H)  that 
are  easily  ionized  in  water.  The  properties  of  hydrated 
Nafion  are  due  to  its  nanophase- segregated  structure  in 
which  hydrophilic  clusters  are  embedded  in  a  hydrophobic 
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matrix.  We  have  characterized  the  nanostructure,  segre¬ 
gation,  position  of  the  ionic  groups,  and  water  distribution 
of  the  hydrated  Nation  membrane  using  non-reactive  FF 
MD  simulations  [4].  It  was  found  that  the  monomeric 
sequence  of  the  Nation  chain  (i.e.  the  block  segregation  of 
hydrophobic  and  ionizible  groups  along  the  linear  polymer 
chain)  affects  the  nanophase- segregation  structure  and 
transport  in  this  material.  We  carried  out  MD  simulations 
of  Nation  using  two  extreme  monomeric  sequences:  the 
first  one  was  very  dispersed  and  the  second  was  very 
blocky  (figure  4). 

The  monomer  sequence  influences  the  water  mobility 
by  affecting  the  membrane  nanosegregation.  We  com¬ 
puted  the  water  diffusion  coefficients  and  found  that  the 
dispersed  (low  blockiness)  membrane  is  —25%  smaller 
than  for  the  blocky  case  (0.46  X  10-5  vs.  0.59  X 
10_5cm2/s  at  300  K).  The  experimental  value, 
0.50  X  10_5cm2/s,  is  within  the  calculated  range.  On  the 
other  hand,  the  vehicular  diffusion  of  hydronium  is  not 
affected  significantly  by  the  monomeric  sequence. 

We  have  also  characterized  the  structure  of  the  water  in 
the  Nafion  membrane  as  a  function  of  the  water  content 
and  found  that  the  hydrophobic  domains  percolate  for 
water  content  above  3-5  water  molecules  per  a  sulfonate 
group.  We  consider  that  this  percolation  of  the  hydrophilic 
domain  is  essential  to  provide  the  pathways  for  proton 
diffusion  responsible  for  the  experimentally  observed 
increase  of  the  proton  conductivity  above  —3  water 
molecules  per  sulfonate  group  in  Nafion. 

Our  atomic  level  characterization  of  the  micro  structure 
of  hydrated  Nafion  and  its  effect  on  the  diffusion  of  water 
and  vehicular  diffusion  of  protons  provides  the  base  for 
proposing  to 

a.  Parameterize  the  ReaxFF  to  fit  QM  calculations  on  Pt, 
hydrated  Nafion,  and  the  interface  between  them  with 
particular  emphasis  on  an  accurate  description  of  the 
proton  transfer  processes.  This  is  a  key  step  toward 
integrating  the  reactivity  and  dynamics  of  the  proton 
migration  into  the  ReaxFF  framework; 
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b.  Characterize  the  pathways  of  proton  diffusion  inside 
the  membrane.  Of  particular  relevance  is  the 
mechanism  of  proton  transfer  through  the  narrow 
channels  connecting  the  hydrophilic  domains  at 
intermediate  water  content.  For  high  water  (>7  H20 
per  S03)  the  proton  transport  depends  little  on  the  S03 
groups  (except  to  provide  the  protons),  however  our 
simulations  suggest  that  at  low  water,  the  sulfonate 
moieties  may  play  an  important  role  in  allowing 
percolation  of  the  proton  networks; 

c.  Using  ReaxFF  compute  the  proton  transfer  rates  as  a 
function  of  the  channel  width  and  sulfonate  distri¬ 
bution.  This  information  is  used  to  develop  a 
mesoscale  model  for  Nafion  membrane; 

d.  Develop  an  electrode  potential  dependent  parameter¬ 
ization  of  the  ReaxFF  to  study  the  electrochemical 
reactions  at  different  applied  potentials; 

e.  Determine  the  mechanism  of  the  proton  transfer  at  the 
interface,  and  the  electrochemical  reactions  at  the 
electrodes;  and 

f.  Determine  water-replacement  polyelectrolytes 
capable  of  high  operational  temperatures  (>  120°C) 
and  high  proton  conductivities  (>0.1  S/cm). 

3.2  Modeling  of  Pt  and  Pt-alloy  electrodes 

Despite  great  effort  to  find  alternative  efficient  catalysts 
for  PEM-FC  applications,  Pt  and  Pt-based  alloys  remain 
the  only  catalysts  capable  of  generating  high  rates  of 
chemical  reactions  at  relatively  low  operating  tempera¬ 
tures  (—  80°C)  of  PEM-FCs.  We  have  used  QM  to  describe 
such  catalysts  in  two  complementary  ways:  (1)  finite 
clusters  containing  —  35  metal  atoms,  and  (2)  infinite  slabs 
having  the  thickness  of  3-5  layers  of  metals,  and  using 
periodic  boundaries. 

The  finite  cluster  calculations  are  more  useful  for 
establishing  the  detailed  sequence  of  reaction  steps 
because  more  accurate  hybrid  levels  of  DFT  are  available. 
In  addition  such  clusters  describe  the  highly  dispersed  Pt- 
nanoparticles  (2-5  nm)  being  used. 
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Figure  4.  Nanostructure  of  hydrated  Nafion.  White  domains  are  mainly  occupied  by  the  Nafion  backbone.  Spheres  represent  sulfur  atoms.  The  surface 
formed  by  dense  dots  is  the  Nafion/water  interface. 
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Figure  5.  Analysis  of  the  occupied  orbitals  for  the  Pt3-CCH2  and  Pt3-CCH  systems. 


Periodic  slabs  are  most  useful  when  comparing  to 
surface  science  experiments  and  for  describing  the 
processes  at  electrodes.  Generally,  we  do  some  calcu¬ 
lations  using  both  approaches  to  allow  crosschecking. 

We  have  carried  out  a  number  of  QM  calculations 
describing  a  Pt-catalyst  using  a  finite  cluster,  increasing 
systematically  the  number  of  cluster-atoms  and  examining 
both  the  electronic  structure  of  the  surface  cluster  and 
reactions,  to  find  the  smallest  suitable  cluster  [7].  Our 
conclusions  from  these  studies  is  that  the  minimum  size 
for  a  cluster  to  accurately  describe  chemisorption  is  the 
three-layer  system  Pt9 10.9  (28  atoms),  which  leads  to 
adsorption  energies,  bond  distances,  and  vibrational 
frequencies  similar  to  the  infinite  Pt(  111)- surf  ace.  To 
describe  surface  reactions  of  organic  molecules  on  the 
catalyst  surface  [e.g.  Pt(lll)]  requires  a  slightly  larger 
cluster  with  35  atoms,  Pt14.13.8-cluster,  allowing  calcu¬ 
lation  of  reactions  on  the  surface  while  minimizing  side 
effects. 

In  addition,  based  on  the  interstitial  electron  model 
(IEM)  of  McAdon  and  Goddard  [22]  as  extended  by  Kua 
and  Goddard  [23],  we  have  developed  an  orbital  model  of 
the  Pt  surface  that  explains  its  reactivity  for  reforming 
various  hydrocarbons  [24].  Figure  5  shows  examples  of 
orbitals  obtained  from  a  Pipek-Mezey  localization  of  the 
occupied  orbitals  for  Pt3-CCH2  and  Pt3-CCH. 


3.2.1  First  principles  study  of  fundamental  reactions 
at  Pt-anode  catalyst.  We  have  already  carried  out  QM 
calculations  on  the  surface  reactivity  of  the  Pt  and  PtRu 
anode  catalysts.  This  QM  uses  a  new  ab  initio  DFT-GGA 
method  (X3FYP)  [6],  which  more  accurately  describes 
such  weak  and  non-bonded  interactions  as  van  der  Waals 
dispersion.  In  order  to  represent  the  surface  of  the  catalyst 
finite  clusters  and  periodic  slabs  will  be  used.  The  overall 
chemical  reaction  at  the  anode  in  a  direct  methane  FC  is: 

CH4  +  2H20  —  C02  +  8H+  +  8e“ . 

We  previously  used  QM  on  small  clusters  to  examine  the 
role  of  various  intermediates.  We  now  extend  these  results 
to  the  35-atom  clusters  (figure  6).  In  addition  to  activating 
the  CH  bond  of  CH4  and  the  OH  bond  of  H20,  what  is 
particularly  important  is  the  final  step  of  forming  the  C02 
product.  We  need  to  avoid  forming  CO,  which  both 
poisons  the  surface  and  reduces  efficiency.  We  are 
particularly  interested  in  the  possible  role  of  oxidized  and 
hydrated  Ru  clusters  being  formed  under  reaction 
conditions.  These  may  play  a  role  in  directing  the  second 
O  to  form  C02  and  in  transporting  the  proton  from  the 
catalyst  to  the  electrolyte. 

To  determine  the  coverage  dependencies  and  the  influence 
of  pre-adsorbed  species,  we  use  periodic  slab-calculations. 


Figure  6.  CH2OH,  CHOH,  COH,  and  CO  intermediates  on  Pt(lll). 
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These  calculations  can  also  include  solvation  effect  through 
the  Poisson -Boltzmann  methods. 

Most  important  here  is  that  the  QM  results  are  used  to 
train  the  ReaxFF  reactive  force  field.  This  allows  us  to 
simulate  very  complex  systems  with  10,000-1,000,000 
atoms,  far  too  large  for  QM,  and  to  examine  the  migration 
of  fuel  to  the  electrode  and  reacting  while  also  treating 
transport  of  the  fuel  through  the  membrane  simultaneous 
with  the  protons  and  water. 

In  addition  to  the  basic  Pt  and  PtRu  systems,  we  plan  to 
explore  how  doping  with  other  elements  might  change  the 
processes  to  minimize  poisons  while  lowering  barriers  to 
the  reactions  and  the  transport  of  protons  and  increasing 
barriers  for  transport  of  fuel.  We  expect  that  this  will  lead 
to  new  strategies  that  can  be  tested  experimentally  to 
develop  improved  anode  catalysts. 

3.2.2  Atomistic  level  modeling  of  cathode  catalytic 
reactions.  The  major  source  of  inefficiency  in  PEM-FC 
remains  to  be  the  significant  overpotentials  (excessive 
reaction  barriers)  in  the  PEM-FC  cathode  reaction  of 
protons  flowing  through  the  membrane  to  convert 
dioxygen  (02)  to  water. 

02  +  4H+  +  4e"  H20, 

(both  for  H2  and  CH3OH  fuel).  This  four-electron 
reduction  process  involves  a  complex  multi-step  sequence 
in  which  there  remains  significant  uncertainty  regarding 
the  reaction  mechanism.  Despite  great  efforts  to  find 
alternative  catalysts  to  promote  a  high  rate  of  oxygen 
reduction,  Pt-based  materials  remain  the  best  for  the 
relatively  low  temperatures  (~  80°C)  at  which  PEM-FCs 
operate.  The  performance  of  PEM-FCs  is  primarily 
limited  by  the  slow  rate  of  the  02  reduction  half  reaction, 
which  is  102  times  slower  than  the  H2  oxidation  half 
reaction.  Therefore,  it  is  very  important  to  understand  the 
reaction  processes  at  the  cathode  in  order  to  find  a 
sequence  that  would  maintain  a  small  barrier  for  each  step. 

We  used  3 5 -atom  clusters  to  determine  the  binding 
energies  and  structures  for  each  possible  intermediate 
(O,  H,  02,  H2,  OH,  OOH,  H20)  involved  in  the  cathode 
reaction  on  the  Pt(lll)  surface  and  some  of  these 
intermediates  on  the  Pt3Ni  and  Pt3Co  alloy  catalysts. 
We  find  that  atomic  oxygen  is  most  strongly  bound  at  the 
p,3-fcc  position,  while  molecular  02  prefers  the  bridge 
site.  OOH  prefers  the  same  geometry  as  02  with  one  O 
covalently  bound  on  top  of  a  Pt  atom  (24kcal/mol). 
Including  zero-point  energy  (ZPE)  a  single  H  atom  prefers 
the  p,3-fcc  site  over  an  on  top  site  by  3.2  kcal/mol,  whereas 
molecular  H2  undergoes  dissociation  to  two  on  top  bound 
H  atoms  wile  adsorbing  on  Pt.  OH  and  water  show 
comparable  binding  structures  (on  top  bound),  but  a 
different  type  of  binding.  The  hydroxyl  radical  is 
covalently  bound  to  one  Pt  atom  (47  kcal/mol),  and 
water  uses  the  remaining  lone  pair  orbital  of  its  oxygen  to 
attach  to  the  surface  atom  (14  kcal/mol). 


In  order  to  determine  the  complete  reaction  pathway, 
we  calculated  all  possible  dissociation  processes  of  the 
various  intermediates  on  the  Pt-cluster.  Using  this 
energetics  we  calculated  heats  of  formation  (A//f)  and 
combined  them  with  the  dissociation  barriers.  Since 
oxygen  and  hydrogen  react  catalytically  at  the  cathode  to 
form  water,  we  determined  possible  reaction  pathways 
starting  with  gas-phase  H2  and  02.  We  distinguished  two 
main  gas-phase  reaction  pathways:  (1)  02-dissociation, 
figure  7(a)  and  (2)  OOH-formation,  figure  7(b). 

Along  the  02-dissociation  pathway  dioxygen  adsorbs 
on  the  surface,  then  dissociates,  and  finally  reacts  with 
chemisorbed  hydrogen  atoms  to  form  water.  The  rate¬ 
determining  step  for  this  mechanism  is  the  Oad  +  Had  — ► 
OHad  reaction  with  a  barrier  of  32  kcal/mol  rather  than  the 
dissociation  of  02,  whose  barrier  is  only  15  kcal/mol. 

Along  the  OOH-formation  pathway,  a  surface  hydro¬ 
gen  first  reacts  with  adsorbed  02  to  form  OOH,  which 
dissociates  to  form  chemisorbed  OH  and  O.  These 
finally  react  with  other  chemisorbed  hydrogen  to  from 
water.  In  this  mechanism  the  OOHad  — ►  OHad  +  Oadc 
dissociation  step  has  the  highest  barrier  of  17  kcal/mol. 
Therefore,  the  OOH-formation  mechanism  seems  favor¬ 
able  for  the  cathode  reaction.  The  results  discussed 
above  were  for  the  gas  phase  processes,  which  we  use  as 
the  starting  point  for  examining  the  cathode  reaction 
mechanisms  under  conditions  suitable  for  fuel  cells.  Fuel 
cells  are  dynamic  systems  in  which  the  catalytic  particles 
at  the  electrode  can  directly  contact  parts  of  the 
membrane  and  surrounding  water  medium.  Thus  we 
must  include  the  effects  of  solvation  and  of  external 
electric  field  (figure  8). 

First,  we  use  DFT  in  combination  with  the  self- 
consistent  reaction  field  (SCRF)  from  Poisson  Boltzmann 
theory  to  describe  the  solvation  effects  for  a  hydrated 
environment  and  determine  the  structural  and  energetic 
changes  due  to  solvation.  For  example,  the  structure  of  a 
layer  of  adsorbed  water  molecules  on  Pt  surface  changes 
drastically  upon  solvation  and  the  binding  energy  to  the 
surface  decreases  significantly.  Polarized  species,  such  as 
chemisorbed  H02  intermediate  show  especially  large 
changes  in  the  structure,  binding,  and  dissociation 
energies  due  to  solvation. 

The  water  environment  also  changes  the  charge 
distribution  of  the  surface  species,  leading  to  both 
increases  and  decreases  in  the  binding  to  the  surface, 
which  changes  the  corresponding  reaction  barriers.  Some 
areas  of  the  catalyst  surface  might  be  hydrophobic  due  to 
nearby  fluorocarbon  regions  of  the  membrane;  these 
regions  may  prefer  02.  Here,  the  reaction  mechanism  may 
resemble  that  in  the  gas-phase.  Other  regions  may  be 
hydrophilic  due  to  nearby  water  channels  and  sulfate 
groups.  Here  solvation  may  drastically  change  the  reaction 
pathways.  Some  reaction  steps  may  occur  in  hydrophobic 
interface  areas  while  the  others  may  prefer  hydrophilic 
areas.  We  must  utilize  both  sets  of  energetics  and 
dissociation  barriers  (gas-phase  and  solution)  to  describe 
fuel  cell  performance. 
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Figure  7.  (a)  Energetics  for  the  02-dissociation 


To  complete  the  picture  of  the  cathode  mechanism,  we 
must  model  the  effect  of  the  electric  field  under  fuel  cell 
operating  conditions.  Thus  for  calculations  in  which  the 
catalyst  is  treated  as  a  slab,  we  add  the  electric  field  to  the 
QM  or  FF  Hamiltonian,  which  has  the  effect  of  modifying 
the  effective  chemical  potential  (Fermi  energy)  as  a 


function  of  the  distance  from  the  metallic  catalyst  surface. 
For  cases  in  which  a  cluster  is  used  to  describe  the  catalysts, 
we  can  treat  the  effect  of  the  field  on  the  reactions  at  the 
catalyst  surface  by  changing  the  net  charge  on  the  catalyst 
(the  gas  phase  calculations  were  carried  for  systems  with 
zero  net  charge).  Calculations  of  adsorbed  water  on  Pt(l  1 1) 
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clusters  with  various  net  charges  show  significant  changes 
in  the  orientation  of  the  water  molecule  (figure  9).  We  will 
carry  out  similar  studies  for  all  intermediates  involved  in 
the  reaction  at  the  conventional  Pt-electrode.  Since  the 
reaction  processes  on  the  Pt-electrode  have  been  well 
studied  experimentally,  we  use  Pt  to  validate  our 
methodologies  for  the  complete  set  of  cathode  reaction 
mechanism,  as  described  above  for  the  gas-phase. 

We  plan  to  perform  the  QM  calculations  on  the  full 
reaction  mechanism  for  the  Pt(lll),  Pt(100),  Pt(110) 
surfaces  and  determine  the  adsorption  energy  of  02  and 
various  intermediates  (OH,  H02)  at  different  terrace  sites 
(hep,  fee,  bridge,  on  top)  and  both  types  of  step  edges. 
Figure  10  shows  02  at  two  such  sites.  We  then  use  QM  to 
examine  the  dissociation  of  the  dioxygen  at  different  sites 
and  the  migration  barriers  between  various  minima.  This 
provides  the  fundamental  elementary  steps  in  the  process 
of  activating  02  and  producing  H20.  Among  the  terrace 
sites  we  find  that  dioxygen  binds  strongest  at  bridge  sites. 
In  contrast,  oxygen  atom  prefers  the  threefold  hollow. 
Clearly  one  must  use  a  model  catalyst  surface  containing 
all  possible  sites  to  describe  properly  the  rates  and  barriers 
for  the  real  catalysts.  We  find  that  such  a  model  requires 
—  300  atoms,  too  large  for  QM  but  quite  straightforward 
for  ReaxFF.  With  ReaxFF  we  can  carry  out  MD  studies  to 
determine  the  effect  of  surface  defects,  step  edges,  kinks, 
impurities,  alloying  elements,  etc.  on  the  barriers  and 
rates. 
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Figure  9.  Structural  changes  due  to  electric  field  effects. 


(a)  (b) 

Figure  10.  Adsorption  of  dioxygen  on  the  fee  site  (a)  and  on  the  [100] 
step  (b). 


An  important  factor  to  include  is  electric  field 
enhancement  [25]  at  step  edges  that  may  amplify  the 
rate  of  the  oxygen  reduction  reaction  (ORR),  the  primary 
factor  limiting  the  PEM-FC  performance.  Using  ReaxFF 
we  can  impose  effect  field  differences  in  different  regions 
by  fixing  the  net  charge  on  the  middle  of  the  cathode 
catalyst  particle  to  be  negative  and  the  middle  of  the  anode 
catalyst  particle  to  be  positive.  Then  the  charge 
equilibration  will  automatically  lead  to  a  potential 
gradient  between  them,  providing  appropriate  field 
enhancement  factors  in  the  observed  rates. 

After  predicting  the  rates  for  the  gas  phase  reactions  we 
can  examine  the  rates  for  the  catalyst  covered  with  water 
and  for  the  case,  in  which  protons  come  directly  from  the 
membrane.  ReaxFF  calculations  for  the  large-scale 
systems  (10,000-1,000,000  atoms)  provide  the  perform¬ 
ance  parameters  as  a  function  of  catalyst  and  membrane 
compositions,  needed  for  estimation  of  the  mesoscale 
parameters  to  characterize  fuel  cell  performance  at 
realistic  conditions. 


3.3  Hydrated  polymer  electrolyte/cathode  interface , 
simulations  of  the  half-cell 

We  believe  that  the  greatest  challenge  to  understanding 
and  optimizing  PEM-FC  performance  is  to  comprehend 
the  processes  taking  place  at  the  interface  between  the 
hydrated  polymer  electrolyte  (Nafion)  and  the  cathode 
(carbon- supported  Pt  nanoparticles).  Since  the  ORR 
occurs  with  protons  transported  from  the  anode  through 
the  membrane  to  the  cathode,  it  is  essential  to  understand 
how  the  protons  are  transferred  across  the  membrane - 
cathode  interface,  particularly,  how  this  is  affected  by  the 
distribution  of  water  at  this  interface.  Although  we  have 
established  that  a  connected  (percolating)  water  structure 
exists  throughout  the  Nafion  membrane,  no  atomistic 
structural  model  has  been  proposed  to  explain  the 
distribution  of  water  and  02  at  the  Pt-Nafion  interface. 

As  the  first  step,  we  propose  to  predict  the  structures  at 
this  interface  using  the  ReaxFF.  In  addition  to  the 
structural  features  of  the  interface,  we  will  determine  the 
dynamics  of  protons,  02  molecules,  and  water  entering 
the  interface  and  reacting.  We  anticipate  that  the  interface 
will  be  heterogeneous  with  the  02  molecule  preferring 
non-polar  hydrophobic  regions,  while  the  proton  (H30+) 
and  H20  prefer  the  hydrophilic  regions.  Thus  the 
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important  chemistry  may  occur  at  the  boundaries  between 
these  regions  and  the  catalyst.  Indeed  determining  the 
distribution  and  dynamics  of  these  important  species  in 
these  regions  will  allow  us  to  examine  how  performance 
might  be  enhanced  by  modifying  the  interface  between  the 
cathode  and  membrane.  This  will  likely  be  most  important 
to  developing  new  PEM-FC  systems  with  improved 
performance. 

To  provide  a  basis  for  answering  such  fundamental 
primary  questions,  we  have  already  carried  out  preliminary 


studies  on  a  model  of  such  an  interface,  including  multiple 
phases  (a  hydrophobic  Nation  backbone,  water,  02,  and 
carbon  support)  in  contact  with  each  other  as  shown  in 
figure  1 1(a).  The  Pt  cathode  nanoparticle  is  modeled  using 
an  87-Pt  atom  cluster,  leading  to  a  diameter  of  —20  A, 
while  the  carbon  support  is  represented  using  graphitic 
layers.  We  constructed  this  interface  by  integrating  it  with 
the  hydrated  Nation  membrane  predicted  previously  [4]. 
Then  this  system  was  equilibrated  using  canonical 
ensemble  MD  simulations  (figure  11(b)).  To  describe  the 


Figure  11.  (a)  Atomic  level  representation  of  the  hydrated  Nafion/cathode  interface,  (b)  Atomistic  model  of  the  cathode  interface.  The  orange  balls  are 

O  atoms  in  02  molecules  and  the  light  blue  balls  are  O  atoms  and  H30+  entities,  (c)  Atomistic  model  used  to  simulate  performance  of  the  overall  PEM 
FC,  including  cathode,  anode,  and  interfaces.  The  orange  balls  are  O  atoms  in  02  molecules  and  the  light  blue  balls  are  O  atoms  in  H30+  entities.  Every  r 
time  steps  an  H2  is  added  between  the  carbon  layers  at  the  top,  every  2r  time  steps  an  02  is  added  between  the  carbon  layers  at  the  bottom,  and  every  r 
time  steps  an  H20  is  removed  above  the  carbon  but  away  from  the  catalysts  at  the  bottom.  The  effective  current  is  one  electron  every  r/2  steps.  The  use  of 
the  ReaxFF  allows  the  reactions  at  the  catalysts  and  the  migration  of  the  protons  to  be  described  with  MD  studies  to  predict  the  overall  fuel  cell 
performance. 
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interaction  of  the  Pt  nanoparticle  with  other  components 
(Nation,  water,  proton  and  oxygen)  in  this  MD  simulation, 
we  optimized  the  FF  to  fit  the  binding  energy  from  DFT 
calculations  at  the  B3LYP/6-31G**  level.  The  FF  for  this 
Pt-cluster  was  optimized  to  reproduce  such  data  as  density, 
structure,  and  heat  capacity.  For  Nation  and  water,  we  used 
the  fluorinated  carbon  FF  [26]  and  the  F3C  water  FF  [27] 
employed  successfully  in  previous  studies  of  the  Nation 
membrane.  This  was  also  used  to  study  a  new  dendrimer 
based  alternative  PEM  material  [28]. 

The  MD  simulation  showed  that  the  hydrated  Nation 
has  a  well-developed  phase- segregated  morphology  with 
percolated  water  distributions  reaching  close  to  the  Pt 
nanoparticle  at  the  cathode.  We  expect  that  such  a 
structure  is  essential  for  proton  transport  to  the  cathode. 
These  simulations  found  that  the  Pt  nanoparticle  is  entirely 
wetted  by  water,  so  that  the  proton  (hydronium)  particle 
can  access  the  entire  Pt  nanoparticle  surface.  On  the  other 
hand,  the  02  molecules  diffuse  reasonably  fast  through  the 
hydrophobic  Nation  backbone,  consistent  with  our 
expectation  that  the  solubility  of  02  is  higher  in  the 
hydrophobic  Nation  backbone  phase  than  in  the  water 
phase. 

Under  constant  current  operating  conditions  water 
molecules  are  generated  while  protons  and  02  molecules 
are  consumed,  all  in  the  presence  of  electric  field.  This  is  a 
highly  non-equilibrium  system,  requiring  use  of  ReaxFF 
to  describe  the  interface  in  such  a  reacting  non¬ 
equilibrium  state.  We  are  currently  carrying  out 
exploratory  calculations  on  this  process.  Since  02 
molecules  and  protons  are  used  up  while  H20  molecules 
are  produced,  we  simulate  the  steady  state  process  by: 

(a)  adding  02  at  a  regular  rate  say  every  AT  steps  between 
the  top  two  layers  of  graphite  but  excluding  the 
footprint  of  the  catalyst  and  turning  off  the 
interactions  of  the  02  with  the  top  graphite  layer; 

(b)  adding  a  proton  to  a  water  at  the  top  of  the  diagram  at 
a  regular  rate  say  every  T  steps; 

(c)  removing  a  water  near  the  graphite  but  away  from  the 
catalyst  at  a  regular  rate  say  every  2 T  steps. 

The  overall  current  is  imposed  by  the  value  of  T 
controlling  the  rates  of  these  processes.  This  process  can 
be  carried  out  for  various  temperatures,  where  we  expect 
that  the  effective  activation  barrier  will  increase  for  higher 
currents. 

We  use  a  similar  strategy  to  describe  the  processes  at  the 
anode -electrolyte  interface  and  to  study  the  properties  of 
this  half-cell  as  a  function  of  current  and  temperatures. 

Combining  both  processes  we  can  describe  the 
performance  of  the  overall  fuel  cell.  For  the  overall  FC 
simulation  (figure  11(c))  the  system  should  be  driven  by 
controlling  the  rates  of  introducing  fuel  (H2  or  CH3OH)  at 
the  anode  and  02  at  the  cathode  and  removing  H20  at  the 
cathode.  These  rates  determine  the  effective  current.  After 
steady  state  is  achieved,  we  calculate  the  net  potential 
difference  between  the  two  electrodes  at  this  fixed  current. 


Varying  the  rates  (current)  changes  the  overall  potential 
leading  to  the  fundamental  performance  figure  of  merit. 
This  will  allow  us  to  determine  how  the  figure  of  merit 
changes  as  the  composition,  configuration,  and  other 
conditions  are  modified. 

4.  Alternative  proton  conducting  membranes 

The  conductivity  of  the  polymeric  membranes  is  water- 
assisted.  To  keep  the  water  content  high,  the  polymer 
electrolyte  must  be  humidified.  This  limits  the  operating 
temperature  of  PEM-FCs  to  90°C  and  adds  further 
complications  to  the  design  of  the  fuel  cell  stack  and 
system.  To  obtain  useful  rates  at  this  low  temperature 
requires  the  very  best  and  most  expensive  platinum 
catalysts,  particularly  for  the  cathode.  The  choice  of  Pt  for 
the  anode  creates  additional  problems  since  CO  can  poison 
a  Pt  catalyst  operating  at  80°C,  even  at  concentrations  of 
only  about  50  parts  per  million.  The  poisoning  could  be 
eliminated  by  increasing  the  operating  temperature  to 
130°C.  However,  the  hydrated  Nafion  loses  too  much 
water  when  it  operates  at  temperatures  above  90°C. 

Another  serious  problem  for  current  polymer  electro¬ 
lytes  is  that  they  are  too  permeable  to  fuel  molecules, 
decreasing  the  efficiency. 

Several  suggestions  have  been  made  for  alternative 
proton  conducting  membranes  to  reduce  the  above 
problems,  including: 

(a)  polymeric  materials  with  heterocyclic  pendant 
groups,  such  as  imidazole,  pyrazole,  or  benzimida¬ 
zole; 

(b)  new  dendrimer-PTFE  copolymers  combining  hydro¬ 
philic  dendrimers  with  hydrophobic  linear  polymers; 

(c)  novel  Brpnsted  acid  based  ionic  liquids,  made  by 
combining  organic  amines  above  their  melting 
temperatures  with  bis(trifluoromethanesulfonyl) 
imide;  and 

(d)  simple  protic  ionic  liquid  mixtures  such  as  a  4:6 
mixture  of  methyl  and  dimethyl  ammonium  nitrate 
that  become  superionic  at  temperature  as  low  as  25°C. 

We  have  already  examined  some  of  these  systems  as 
described  below. 

4.1  Imidazole-based  water-free  proton  conducting 
membranes 

As  mentioned  above,  the  very  high  protonic  conductivity 
of  the  hydrated  perfluorosulfonic  polymers  is  due  to 
percolating  water  channels  surrounded  by  sulfonic  groups 
to  protonate  the  water.  The  proton  conductivity  for  such 
systems  is  quite  high,  ~10-1S/cm  for  Nafion  at  85°C. 
However,  substituting  water  with  such  heterocycles,  as 
imidazole,  pyrazole,  and  benzimidazole  has  led  to  some 
cases  with  comparable  conductivity  [29]. 

Imidazole  has  a  high  boiling  point  and  good  proton 
transfer  capability,  suggesting  it  as  a  candidate  to  replace 
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water  as  a  proton  transfer  molecule  in  Nation  membranes 
[30] .  Indeed  imidazole  provides  a  very  good  conductivity, 
~10_1S/cm  at  160-180°C.  Unfortunately  it  was  found 
that  imidazole  chemisorbs  on  the  Pt  electrode,  poisoning 
the  catalyst. 

Recently,  we  showed  that  fluorinated  imidazole  in 
Nation  [20]  are  an  excellent  candidate  for  a  high 
temperature  PEM,  combining  both  high  boiling  point 
and  conductivity  required  for  the  electrolyte  while  not 
chemisorbing  on  the  electrode.  We  used  QM  to  show  that 
fluorinating  imidazole  dramatically  decreases  the  binding 
energy  of  imidazole  to  Pt  (from  21.1  to  1.3kcal/mol  for 
trifluoroimidazole),  to  a  value  much  lower  than  the 
binding  energy  for  water  on  Pt  (12.5  kcal/mol). 

Moreover,  we  calculated  the  proton  diffusion  in  the 
trifluoroimidazole/Nafion  membrane  to  be  only  40% 
lower  than  that  of  imidazole/Nafion  at  the  same 
temperature  and  concentration.  This  comparison  of  the 
proton  diffusivities  was  done  by  a  combined  QM/MD 
approach  we  developed  for  such  calculations.  In  this 
approach,  we  compute  the  vehicular  diffusion  component 
from  the  classical  MD  simulations  using  non-reactive  FF 
with  explicit  H30+  entities.  Then  we  calculate  the 
hopping  diffusion  term  using  the  configurations  from  the 
MD  and  evaluating  the  proton  hop  rates  using  transition 
state  theory  (TST)  with  quantum  tunneling  correction 
based  on  energy  curves  for  proton  transfer  completely 
computed  from  ab  initio  QM.  Thus  trifluoroimidazole  is  a 
promising  candidate  for  water  replacement. 

The  method  used  above  to  estimate  proton  diffusion 
rates  in  the  bulk  electrolyte  is  simple  and  efficient  way,  but 
it  is  not  suitable  for  describing  details  of  proton  transfer  at 
the  electrolyte -catalyst  interface.  We  intend  to  use 
ReaxFF  MD  methodology,  allowing  explicit  proton 
hopping  to  obtain  the  rates. 

4.2  Dendrion  polymeric  materials  combining 
hydrophilic  dendrimer  and  hydrophobic  PTFE 

Current  state-of-the-art  PEM-FC  technologies  are  based 
on  the  hydrated  perfluoro sulfonic  polymers  such  as 
Nafion,  which  has  a  high  proton  conductivity  as  well  as 
excellent  chemical  and  thermal  stability.  These  favorable 
properties  result  from  the  nanophase- segregated  structure 
consisting  of  two  phases:  one  is  the  hydrophobic  phase  of 
perfluorinated  polymer  backbone  and  the  other  is  the 
hydrophilic  water  phase  associated  with  sulfonic  groups 
belonging  to  polymer.  In  previous  studies  [4],  we  found 
that  the  nanophase- segregated  structure  can  be  controlled 
by  specifying  the  molecular  architecture  (monomeric 
sequence  and  cross  linking),  which  in  turn  determines  the 
transport  characteristics  of  PEMs. 

Recently  we  proposed  the  Dendrion  copolymer  in 
which  a  water-soluble  poly  aryl  ethereal  dendrimer 
(Frechet  type  dendrimer  [31])  is  covalently  bound  with  a 
hydrophobic  linear  PTFE.  From  such  molecular  constitu¬ 
ents  we  expected  that  this  copolymer  would  form  a 
nanophase-segregated  structure  with  a  specific  amount  of 


water  available  for  proton  transport.  We  also  expected  that 
the  polar  groups  were  localized  at  the  surface  of  the 
dendrimer  through  covalent  connectivity,  leading  to 
formation  of  a  well-defined  spatial  distribution  of  polar 
group  in  the  membrane.  Since  the  transport  properties  of 
the  membrane  are  coupled  with  its  nanostructure,  we  may 
systematically  control  the  transport  properties  by 
managing  water  distribution  through  nanophase  segre¬ 
gation  in  the  membrane.  Many  factors  affect  the 
nanophase  segregation  (dendrimer  size,  type  of  acid 
groups,  equivalent  weight  of  copolymer),  but  our  initial 
studies  focused  on  the  water  content  and  distribution  as 
a  function  of  temperature,  since  this  is  critical  to  good 
performance. 

We  carried  out  MD  simulations  to  determine  the 
structure  and  properties  of  the  Dendrion  hydrated 
copolymer  membrane  consisting  of  the  second-generation 
Frechet  type  dendrimer  combined  with  linear  PTFE 
copolymer  (1100  of  equivalent  weight)  to  assess  the  effect 
of  water  content  on  the  nanophase- segregated  structure 
and  transport  properties. 

We  prepared  two  model  systems  each  of  the  same 
dendrimer-PTFE  copolymer  chains  but  with  different  water 
contents  (5  and  10wt%).  Then  we  carried  out  MD 
simulations  under  the  conditions  previously  used  for  Nafion: 
1  atm  and  353.15  K.  We  found  that  nanophase  segregation, 
which  appeared  during  the  pressure -temperature  annealing 
procedure,  used  to  determine  the  equilibrium  structure  of 
such  amorphous  polymers  (figure  12).  We  also  examined  the 
structure  and  the  transport  properties  from  a  5  ns  equilibrium 
NPT  MD  simulation. 

Analyzing  the  water  structure  in  the  Dendrion,  we 
found  that  the  water  phase  formed  bulk-water-like 
connected  channels  more  efficiently,  than  in  the  Nafion 
membrane.  We  expect  that  this  bulk-water-like  structure 
will  facilitate  proton  hopping  as  easily  as  for  the  bulk 
water. 

We  calculated  a  water  diffusion  coefficient  of 
0.37  X  10_5cm2/s  (10wt%  water  content)  and  0.21  X 
10_5cm2/s  (5wt%  water  content),  which  is  4-7  times 
smaller  than  in  Nafion  membrane  (—1.5  X  10_5cm2/s). 
This  suggests  that  the  Dendrion  may  have  a  smaller 
electro-osmotic  drag  coefficient  than  Nafion. 


Figure  12.  Predicted  nanophase  segregated  structure  of  the  dendrion 
hydrated  Frechet  dendrimer-PTFE  copolymer  membrane. 
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Figure  13.  Proton  transfer  barriers  as  a  function  of  R(H-  •  O)  distance 
(in  A)  in  the  methylethylammonium  nitrate  ionic  liquid. 

These  structural  and  dynamical  features  of  the  dendrion 
membrane  make  it  a  promising  candidate  for  FC 
applications  [28]. 

4.3  Ionic  liquids 

Just  recently  it  was  recognized  that  ionic  liquids  can  serve 
as  proton  conducting  non-aqueous  electrolytes  in  FCs. 
Novel  Brpnsted  acid  based  ionic  liquids  made  by 
combining  imidazole  with  bis(trifluoromethanesulfonyl) 
imide  (HTFSI)  are  promising  candidate  materials  for  fuel 
cell  applications  [32].  They  are  electroactive  for  H2 
oxidation  and  02  reduction  at  the  Pt  electrode  under  non¬ 
humidifying  conditions  at  moderate  temperatures 
(~130°C).  Other  very  promising  materials  are  simple 
protic  ionic  liquid  mixtures  (such  as  a  4:6  mixture  of 
methyl  and  dimethyl  ammonium  nitrate)  that  become 
superionic  at  25°C.  At  25°C  the  conductivity  is  similar  to 
Nation  at  80°C  (0.15  vs.  0.1  S/cm),  but  at  100°C  the 
conductivity  of  the  ionic  liquid  system  is  of  0.47  S/cm,  4 
times  larger  than  for  Nation. 

In  order  to  provide  a  basis  for  understanding  the 
superpro tonic  conductivity  of  ionic  liquids,  we  determined 
the  QM  barrier  for  proton  transfer  in  the  methyl -ethyl 
ammonium  nitrate  system  (figure  13).  For  hydrogen  bonds 
shorter  than  2.5  A  there  is  essentially  no  barrier  for  proton 
transfer,  which  is  likely  an  important  factor  in  explaining 
the  superionic  nature  of  this  ionic  liquid. 

For  these  systems  to  be  used  in  a  commercial  FC,  we  want 
to  anchor  the  anion  of  the  ionic  liquid  to  an  immobile  phase, 


preferably  a  high  temperature  stable  polymer.  A  variety  of 
ionic  liquids  (such  as  in  figure  14)  will  be  screened  and  then 
the  most  promising  candidates  will  be  selected  for  further 
study.  We  consider  that  such  anhydrous  protonic  materials 
offer  excellent  opportunities  for  improved  FC  materials.  It  is 
important  to  determine  the  fundamental  physics  and 
chemistry  of  the  proton  conducting  phenomena  in  this  vast 
under  explored  class  of  materials.  This  will  provide  critical 
data  needed  to  design  protonic  materials  that  can  operate  at 
moderate  temperatures  (100-200°C),  to  attain  significant 
gains  in  efficiency  (conductivity),  low  toxicity  and  high 
thermal  stability,  and  compatibility  with  the  electrode 
catalyst.  We  will  implement  multi-scale  tools  for  quantum 
hopping  molecular  dynamics  (QHMD)  methods  to  simulate 
both  quantum  and  classical  events  in  proton  transfer 
reactions,  including  new  first  principles-based  FFs  capable 
of  proton  transfer  events. 

5.  Electrode  alloy  catalyst  materials  for  PEM-FCs 

5.1  Atomic  scale  design  of  advanced  alloy  anode 
catalysts 

Although  platinum  is  currently  the  best  anode  electro¬ 
catalyst  for  PEM-FCs,  it  suffers  from  CO  poisoning. 
Therefore,  it  is  highly  desirable  to  develop  new  anode 
catalysts  with  higher  CO  tolerance.  Ruthenium  (Ru) 
alloyed  Pt  with  the  atomic  ratio  of  1:1  is  known  to  be 
effective  in  CO  removal  and  to  significantly  promote  FC 
performance.  In  addition,  the  Pt50Ru5o  alloy  on  a  PEM  is 
currently  the  best  anode  catalyst  for  direct  methanol  fuel 
cell  (DMFC)  [33-35]. 

The  mechanism  widely  accepted  for  the  CO  oxidation 
[11,36]  on  PtRu  is: 

Ru  +  H20  — ►  Ru-OH  +  Hadd; 

Pt— CO  +  Ru— OH  — *  Pt  T  Ru  +  C02  -1-  Hadd, 

where  the  rate-determining  step  is  the  formation  of  Ru-OH 
species  [37].  The  relatively  easy  formation  of  hydroxides 
(RuOH)  or  hydrous  oxides  (RuO^Hy)  of  oxyphilic  Ru 
probably  plays  a  key  role  in  preventing  the  CO  poisoning  of 
Pt.  Ru  is  well  known  [11]  for  its  water-dissociating 
capability  and  rich  oxygen  chemistry.  Thus  chemisorbed 
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Figure  14.  Examples  of  proton  acceptor  for  a  Brpnsted  acid  bis(trifluoromethanesulfonyl)  for  non-aqueous  proton  transfer  membrane  media. 


Fuel  cell  catalysts  and  membranes 


265 


oxygen,  surface  oxide,  subsurface  oxygen,  and  buried 
oxides  may  coexist  in  the  near-surface  region  [38,39] .  It  has 
been  reported  that  Pt  remains  metallic  under  operating 
conditions,  while  Ru  forms  oxides,  hydroxides,  or 
amorphous  hydrous  oxides  during  the  electrocycle  [40,41]. 

However,  there  remain  controversies  on  the  exact  nature 
of  the  promotion  mechanism,  both  on  the  nature  of  active 
Ru  species  [metallic  ruthenium  (Ru),  ruthenium  oxides 
(Ru02  and  Ru03),  hydroxides  (RuOH)  or  hydrous  oxides 
(RuOJTy)],  and  also  on  how  the  OH  species  bound  to  Ru 
more  strongly  than  to  Pt  can  spillover  to  a  neighboring  Pt 
sites  and  combine  with  CO  there.  Possibly  it  picks  up  a 
hydrogen  from  a  water  molecule  on  a  neighboring  Pt 
atom,  effectively  a  transfer  of  OH  from  Ru  to  Pt  [42]. 
Therefore,  it  is  important  for  designing  the  most  active 
Pt-Ru  electrocatalysts  to  have 

•  the  optimum  Pt-Ru  surface  concentration, 

•  the  optimum  Pt-Ru  distribution  and  structure, 

•  the  optimum  water  content  around  the  active  sites. 

The  surface  composition  and  structure  of  the  alloy 
catalyst  affect  the  activity  and  selectivity.  This  depends 
drastically  on  such  experimental  conditions  as  the 
synthesis  temperature  and  subsequent  processing  con¬ 
ditions  that  make  complicate  to  extract  an  atomistic 
interpretation  of  the  experimental  activity.  We  use  theory 
and  computation  to  interpret  experiments  by  determining 
the  atomic  level  reaction  mechanisms  and  energetics  for 
various  plausible  surface  structures  and  compositions, 
including  the  boundary  between  the  Pt-Ru  solid  solution 
phase  and  the  hydrous  ruthenium  oxide  phase. 

We  investigate  the  oxygen  and  water  chemistry  of 
hydrous  Ru  and  RuOx  surfaces  (step  1 :  Ru/RuOx/RuOxHy), 
and  then  mix  them  with  Pt  (step  2:  Pt-Ru/RuOx/RuOxHv) 
in  order  to  pin  down  the  active  sites  for  oxygen/water 
activation  and  for  CO  oxidation  to  C02.  This  allows  us  to 
deduce  which  arrangement  of  Pt,  Ru,  and  their 
oxides/hydroxides/hydrous  oxides,  and  possibly  the  other 
alloying  elements,  can  form  appropriate  active  sites  to 
reduce  the  energy  barrier  of  the  critical  steps.  We  use  the 
first-principles  calculations  to  elucidate: 

(a)  surface  compositions  and  structures  of  metal/metal 
oxide  alloys,  including  surface  defects,  for  various 
environments  (oxygen,  water,  etc.;  surface/subsurface 
oxide,  hydroxide,  or  hydrous  oxide  formation); 

(b)  each  step  from  activation  of  the  fuel  (H2)  and  the 
poison  (CO)  and  how  this  depends  on  surface  defects 
and  alloying; 

(c)  how  each  of  these  steps  is  affected  by  solvation 
and  electric  fields. 

All  previous  theoretical  studies  have  been  limited  to 
ideal  metal  surfaces  or  ideal  surfaces  of  Pt-Ru  solid 
solution  phases  [43-48].  We  intend  to  simulate  the 
evolution  of  “realistic”  Pt-Ru  alloy  surfaces  at  operating 
conditions,  including  surface  hydration  and  oxidation. 


It  is  important  to  study  the  surface  segregation  of  Pt,  how  the 
segregation  changes  under  a  variety  of  oxidation  conditions, 
the  effect  of  lattice  strain  (compressive  or  expansive)  as  a 
result  of  the  alloying  process,  and  the  electronic  coupling  of 
the  Pt-skin  layer  with  the  alloy  substrate.  This  should  provide 
insight  on  the  fundamental  origin  of  the  reactivity  difference 
between  the  alloys  and  pure  Pt. 

Since  CO  poisoning  represents  one  of  the  major  problems 
for  the  catalysis  on  Pt-DMFC  anodes,  we  plan  to  study  the 
effect  of  alloying  (PtRu,  Ptlr,  PtOs,  PtMo)  on  the  MeOH 
chemistry.  Our  goal  is  to  develop  insight  into  the  effect  of  the 
electronic  structure  on  specific  bond  breaking  events,  which 
would  be  useful  for  guiding  the  catalyst  discovery  to  more 
promising  areas  of  the  complicated  alloy  phase  diagrams. 
Similarly,  our  DFT  calculations  can  provide  detailed 
information  on  the  adsorbate/alloy-induced  change  in  the 
electronic  structure  of  the  catalytic  surface,  the  structure  of 
surface  reactive  intermediates,  and  diffusion  barriers. 

Along  with  our  studies  of  the  MeOH  chemistry  on  Pt 
and  its  alloys,  we  probe  the  effect  of  water  on  that 
chemistry  as  a  function  of  the  electronic  structure  of  the 
substrate.  Water  is  present  at  the  anode  of  DMFCs  under 
realistic  operating  conditions  and  this  part  of  our 
investigations  will  address  key  mechanistic  issues  of  the 
anode  reactivity.  H20  is  expected  to  alter  the  interaction  of 
MeOH  and  its  derivate  reaction  intermediates  with  the 
catalytic  surface  in  a  substantial  way.  Besides  changing 
the  potential  energy  characterizing  the  MeOH/Pt  inter¬ 
action,  new  reaction  paths  would  become  available  to 
these  intermediates.  One  of  the  most  interesting  parts  of 
this  research  is  to  determine  the  interaction  of  CO  and 
H20  at  the  anode.  As  MeOH  represents  the  first  member 
of  hydrocarbons,  which  can  be  reformed  through  the 
aqueous  phase  hydrocarbons  reforming  (ACR)  process  to 
hydrogen,  our  studies  of  the  MeOH  reactivity  on  the  anode 
catalyst  would  directly  benefit  our  studies  of  the  ACR 
process.  Pt  represents  one  of  the  most  promising  catalysts 
for  the  anode  of  DMFCs  and  for  the  ACR  process. 
Therefore,  reactivity  of  oxygenates  on  Pt  and  Pt-based 
alloys  would  have  an  impact  on  both  subjects. 

We  also  plan  to  study  the  MeOH  decomposition  on 
stepped  Pt  surfaces,  which  is  critical  for  addressing 
particle  size  effects  relevant  to  the  anode  catalysis  of 
DMFCs.  In  particular,  the  mechanism  of  MeOH 
decomposition  will  be  studied  on  the  Pt(2 11)  surface, 
presenting  a  more  reactive  step  edge,  which  could  serve  as 
a  model  for  defect  sites  on  the  catalyst. 

5.2  Atomic  scale  design  of  advanced  alloy  cathode 
catalysts 

The  PEM-FC  performance  is  limited  most  by  the  ORR 
kinetics  and  second  by  transport  losses  due  to  proton 
conduction  in  the  cathode  catalyst  layer  and  oxygen 
transport  to  the  cathode  catalyst.  Therefore,  optimization 
of  the  electrocatalytic  processes  at  the  cathode  could 
significantly  improve  fuel  cell  performance.  Some  results 
showing  activity  gains  have  been  reported  in  the  literature 
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on  binary  alloy  catalysts,  such  as  PtNi,  PtCo,  and  PtCr  [49] 
as  well  as  ternary  alloy  catalysts  [50].  However,  these  have 
not  been  implemented  widely  in  FC  stacks  yet. 

To  make  progress  in  developing  improved  catalysts,  it  is 
important  to  understand  the  fundamental  processes  for 
these  catalysts.  We  studied  Pt3Ni  and  Pt3Co,  since  both 
provide  improved  performance  relative  to  pure  Pt.  We  find 
that  the  bulk  alloy  has  only  Pt  neighbors  to  each  Ni  or  Co, 
which  for  the  surface  would  lead  to  1/4  of  the  surface 
atoms  being  Ni  or  Co.  However,  QM  calculations  for  the 
two-dimensional  periodic  slabs  shows  that  the  most  stable 
surface  has  no  Ni  or  Co  on  the  surface.  Apparently  the 
preference  for  Ni  or  Co  to  have  12-Pt  neighbors  is  strong 
enough  to  drive  what  would  have  been  surface  Ni  or  Co  to 
the  second  layer.  The  result  is  100%  Pt  in  the  top  layer, 
50%  in  the  second  layer  and  then  75%  in  the  third  and 
subsequent  layers.  This  has  also  been  observed  exper¬ 
imentally.  Although  there  are  no  Ni  or  Co  in  the  top  layer, 
the  chemical  activity  of  the  surface  Pt  is  dramatically 
affected  by  where  the  second  layer  Ni  or  Co  are  located. 
This  is  shown  in  figure  15,  where  one  can  see  that  the 
energy  surface  for  movement  of  chemisorbed  O  and  H  is 
quite  corrugated  [5],  with  H  surface  coverage  much  higher 
barriers  than  for  pure  Pt.  Oxygen  has  a  hopping  barrier  of 
~  2  eV,  compared  to  ~  0.6  eV  for  the  pure  Pt( 111)  surface. 
Thus  on  the  alloy  surface  oxygen  is  strongly  localized. 
Similarly,  hydrogen  is  known  to  be  very  mobile  on  the  Pt 
surface,  with  diffusion  barriers  of  only  —0.1  eV.  However, 
for  the  Pt/Ni  alloy  we  found  that  the  hydrogen  atoms  were 
much  less  mobile  leading  to  “channels”  between  and  on 
circles  around  the  Ni  atoms  (figure  15).  The  nearly  fixed 
oxygen  location  and  the  constrained  hydrogen  mobility 
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may  increase  the  probability  that  both  atoms  come  close 
enough  to  react  with  each  other.  This  example  shows  how 
alloying  significantly  changes  the  system  character  and 
has  a  dramatic  effect  on  the  barriers  for  the  various 
reaction  steps  discussed  in  figure  7.  In  particular,  we 
expect  that  the  barrier  for  02  dissociation  will  decrease 
because  of  charge  transfer  from  second  layer  Ni  or  Co  to 
the  surface  Pt,  making  formation  of  chemisorbed  02  more 
favorable. 

We  have  also  studied  finite  clusters  for  the  Pt3Ni  alloy 
and  found  these  results  to  be  quite  different  than  in  the 
infinite  slab.  Ni  atoms  are  found  at  the  surface  with  the 
same  distributions  as  if  they  were  cut  from  the  infinite 
crystal.  Therefore,  clusters  are  not  large  enough  for  the  top 
layer  to  be  pure  Pt.  For  clusters  less  than  4  nm  on  a  side, 
we  expect  the  surface  chemistry  to  be  quite  different  than 
for  the  bulk. 

The  QM  calculations  described  above  are  essential  to 
determine  accurate  barriers  and  pre-exponential  factors. 
However,  for  simulating  complex  processes  on  realistic 
surfaces  as  a  function  of  temperature  and  the  dependence 
on  the  nature  of  the  catalyst -electrolyte  interface,  we 
must  carry  out  MD  simulations  for  at  least  ns  on  systems 
with  many  thousands  of  atoms.  We  extend  ReaxFF  to 
describe  the  various  cathode  catalytic  processes  on  Pt 
alloy  systems. 

6.  Summary 

A  great  deal  of  progress  has  been  made  over  the  last  decade 
in  developing  fuel  cells  with  increased  efficiency  and 
decreased  cost,  but  the  current  generation  remains  well 
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Figure  15.  (a)  Energy  surface  (eV)  for  migration  of  chemisorbed  H  atom  on  the  Pt3Ni(lll)  alloy  surface;  (b)  same  for  O  atom. 
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short  of  the  needs  of  society.  We  propose  to  address  such 
problems  using  recently  developed  rational  approaches 
based  on  first  principles  predictions.  This  uses  the  new  first- 
principles  based  reactive  ReaxFF  force  field  to  describe 
chemical  processes  (ion,  proton  and  water  transfer,  oxygen 
reduction,  reactions  of  intermediates,  phase  transform¬ 
ations,  etc.)  at  a  scale  suitable  to  address  the  issues  of  anode 
catalysis,  cathode  catalysis,  proton-conducting  polymer 
electrolytes,  and  electrode/electrolyte  interfaces.  We  focus 
on  overcoming  some  key  bottlenecks  in  this  process  and  are 
now  at  the  point  where  building  these  techniques  into  an 
integrated  computational  facility  and  applying  them  to  the 
fundamental  processes  at  the  anode  and  cathode  catalysts, 
in  the  proton-conducting  electrolyte  connecting  them,  in 
the  interfaces  between  these  elements,  and  in  the  fuel 
processing  would  enable  the  conception,  synthesis, 
fabrication,  development,  and  understanding  of  advanced 
materials  and  structures  for  PEM-FCs. 
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